Quantum motion effects in an ultracold-atom Mach-Zehnder interferometer 
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We study the effect of quantum motion in a Madi-Zefmder interferometer where uhracold, two- 
level atoms cross a 7r/2-7r-7r/2 configuration of separated, laser illuminated regions. Explicit and 
exact expressions are obtained for transmission amplitudes of monochromatic, incident atomic waves 
using recurrence relations which take into account all possible paths; the direct ones usually con- 
sidered in the simple semiclassical treatment, but including quantum motion corrections, and the 
paths in which the atoms are repeatedly reflected at the flelds. 
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I. INTRODUCTION 

The fringes of an interferometer are sensitive to differ- 
ential phases of the arms caused by unequal fields along 
the interfering paths. This makes interferometers use- 
ful for metrology and fundamental studies. In particu- 
lar, atom interferometers offer, because of the internal 
structure of the atom, richer interactions, greater and 
simpler control than the ones based on light, electrons 
or neutrons [i|. They are used for the precise deter- 
mination of frequencies and times in atomic clocks, as 
well as many other applications to measure with unprece- 
dented accuracy the gravity field '^ and gravity gradients 
[3| , rotations [j, Q , fundamental constants @ , accelera- 
tions [3], or relativistic effects Q. Indeed, the accuracy 
level is currently so high that the theoretical treatments 
of the global performance of the interferometer [3, [i3| 
or the individual constituents (beam splitters, mirrors) 
[ir, ^JZ\ need to be refined with respect to the simple, 
original modelings. A further reason is the use of ultra- 
cold atoms [a, llJ, HJ] to minimize velocity broadening 
and to increase coherence lengths and flight times; they 
also make possible the spatial separation of the arms by 
atomic recoil [15l [iGj , as in Sagnac interferometry, where 
slower atoms increase arm separation, the area enclosed, 
and thus the sensitivity achieved. In addition, the use of 
condensates and other ultracold-matter phases (such as 
the Tonks- Girardeau gas) in internal-state interferome- 
try is currently being explored [ij, |T^ . Since the atomic 
velocities may be nowadays several orders of magnitude 
smaller than in early beam experiments [19], these devel- 
opments raise the following question: Is there any fun- 
damental or practical lower bound for the velocities in 
interferometry? [20| • To answer it we need to go beyond 
the approximation in which the center of mass motion 
along the interferometer arms is treated classically. 
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FIG. 1: Atom interferometry with a 7r/2 — tt — 7r/2 pulse 
configuration. The left edges of the three pulses are located 
at a;i = 0, X2 = 1/ + //2 and xa = 2L + 3Z/2. The beams are 
not drawn to scale and their width is greatly exaggerated. 



Much work on that line by Borde and coworkers has 
emphasized a wave packet approach 9, 11]. We shall ex- 
plore here a complementary stationary method, extend- 
ing some previous results on recurrence relations which 
were applied to Ramsey interferometry in a waveguide 
|2]| . The analysis of stationary solutions leads to use- 
ful insight, and quite frequently provides sufficient in- 
formation, as the history of scattering theory demon- 
strates (consider, e.g., the cavalier but straightforward 
derivation of cross sections from stationary waves versus 
the more rigorous and cumbersome, but finally equiv- 
alent, wave packet derivation). Of course, wave pack- 
ets can be constructed afterwards by linear superposi- 
tion for examining transients and specific space-time pro- 
cesses. Among the interferometers with spatially sepa- 
rated paths we shall focus on the simplest configuration, 
a Mach-Zehnder interferometer, first implemented in the 
time domain by Kasevich and Chu [2]. We shall deal 
here with the version in which the laser beams are sep- 
arated in space [51. It consists on a first 7r/2 laser beam 
acting as an atom-beam splitter, followed by a mirror (tt 
beam) and finally a second, recombining 7r/2 beam, see 
Fig. [TJ Our main tool in this investigation is the im- 
plementation of exact relations for the final transmission 
amplitude in the excited state. They may be cast as "re- 
currence relations" in terms of the scattering amplitudes 



for each laser field [2l|, |2^, which allows us to classify 
and calculate all possible paths by the number of reflec- 
tions ^: the dominant or "direct" ones (without reflec- 
tions), associated with the usual semiclassical ordering 
of events but affected by quantum corrections, and also 
those paths in which the particle is reflected in several 
field regions. The extreme low-velocity regime in which 
these later "multiple scattering" paths become significant 
distorts severely the interference pattern and thus sets a 
fundamental lower limit to the atomic velocities for in- 
terferometry with fields separated in the space domain 
[201 l2l|; for intermediate velocities, just above the mul- 
tiple scattering regime, direct paths dominate, but the 
semiclassical expressions are not yet quite accurate and 
need correction. Therefore, an understanding of the var- 
ious effects and scales involved is useful. To simplify the 
analysis and isolate quantum motion effects from other 
phenomena we shall ignore in this paper any external 
fields different from the laser fields. Other simplifying as- 
sumptions are the consideration of flat and sharp (square) 
laser sheets, fully coherent processes (i.e., we neglect ex- 
cited state relaxation), and semiclassical atom-laser in- 
teraction. Some of these approximations are discussed in 
the final section. 



II. NOTATION AND HAMILTONIAN 



Atom field interaction in 3D 



counter propagating lasers which induce a two-photon 
Raman transition and a large (optical) recoil and arm 
separation, so that the parameters are effective ones, af- 
ter adiabatic elimination of a non-resonant upper state, 
see e.g. [2j, [2J]. In a field adapted interaction-picture 
defined by Hq = hujL\e){e\, and applying the rotating- 
wave approximation (RWA) , the time dependence of the 
Hamiltonian is removed, 
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e^lkLV^H^)] _^ jj^^ 



(2) 

where Aq = lul — ^ge is the detuning between the laser 
frequency and the internal transition. 



B. ID effective equation in s-direction 

To solve the stationary Shrodinger equation 



Ek\4>k{x,y,z)) for an energy 
Ek = h^P/{2m), with P = kl + k"^ + k^, we use the 
ansatz 



\M^, y, z)) = g,{x)e'''yye'''''\g) + e,(x)e*«='^e''=^^|e), 

(3) 
which describes the momentum transfer in y-direction 
when the atom is excited, Qy = ky + k^, and conservation 
of momentum (free evolution) in z-direction. Inserting 
this ansatz into the Shrodinger equation gives an effective 
equation in x direction. 



We consider a setup where a two-level atom, with 
an internal (hyperfine) transition frequency LOge between 
levels \g) and |e), moves with an initial wavenumber 
k = {kx,ky,kz) {kx ^ ky,kz) and is illuminated in three 
x-localized regions by a classical electric field E{x,t) = 
Eq(x) cos [ujLt — k^y + 4>{x)) traveling in y-direction, see 
Fig. [TJ The full 3D Hamiltonian describing this system 
in the Shrodinger picture is 

H = — + nujge\e){e\ 

+ hQ,{x){a++ (J -) cos [ujLt- kLy + 4>{x)], (1) 

where ct+ — \e){g\, ct_ = |ff)(e|, the x-dependent Rabi 
frequency ri{x) is assumed to be constant inside the field 
regions, Q.{x) = f} for x € [0, 1/2], x€\L + 1/2, L + 3^/2] 
and x € [2L -\-Zl/2,2L-\- 21] and zero otherwise, and the 
laser phase (t){x) is constant within each of the illumi- 
nated regions with values 4>m n = 1,2,3. In practice the 
traveling wave is an effective one corresponding to two 
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and A is the effective detuning, 

A = Ao-Afe,, 



Afei< 
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which includes the ordinary detuning Aq and the "kinetic 
detuning" Afcj„ with a photon recoil term and a Doppler 
term. From now on we shall deal with the ID effective 
equation ([4]) only. 



III. SEMICLASSICAL REGIME 



^ In this paper the term "reflection" refers to a change in the sign 
of the momentum component in the longitudinal x-direction. Do 
not confuse this with the recoil taking place at the second (mir- 
ror) laser, in which the excited or ground state are interchanged 
but the momentum component in x-direction does not change 
sign. 



In the simplest treatment, valid for fast enough par- 
ticles, Ex ^ fi^, kx ^ Qx = \/kx + 2mA/h, the in- 
ternal and longitudinal degrees of freedom are decou- 
pled, and the x-component of the center of mass is as- 
sumed to follow the classical trajectory x(t) = Vxt, where 
Vx = hkx/m. We shall, in other words, treat the interfer- 
ometer with fields separated in space as an interferometer 



for fields separated in time. Then, the internal states will 
evolve with H„ , 



Hn = -hA\e){e\ 



nn 



(e-^"a+ 



H.c) 



(8) 



in the nth laser field, and with the bare Hamiltonian 



HB = -nA\e){e\ 



(9) 



in the non-interacting regions. The corresponding time 
evolution operators are e"*^"*/'' in the nth laser and 
^-iHst/h jj^ ^jjg £j,gg evolution regions. From Fig. [l]it can 
be seen that there are four possible semiclassical paths 
which lead to an excited atom from an atom which is 
initially in the ground state \g). If we set r = l/vx and 
T — L/vx, the amplitudes of these four paths are given 

by 
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'iH3r/{2h)-iHBT/n\ 



(e|e — 'v--e --^ ' ■'\g){g\e-'''^^^^e~'''^' '''\e){e\e-''''^^^'^^\g), 



-iH3T/{2h) -iHsT/hi 



(ele —'-■-£ -^^'-\e){e\e-'"'^^^e-'"^'^^\g){g\e-'"'^^'-^''>\g), 



'iH3T/{2h)^iHBT/hi 



(ele -^''^-'e — ■%)(g|e-'«^-/'^e-^«-^/''|g)(g|e-^«^^/^^'^)|g). 



(10) 
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For a wave-packet, the momentum recoil will lead to a 
separation of these paths in the y-direction. If this sepa- 
ration is larger than the transversal position spread of the 
packet Ay, hk^T / {2m) ^ Ay, and the detector's resolu- 
tion is better than hkLT/(2m), the interference between 
the outer paths will be suppressed and only the inter- 
ference between Al'^' and A^'^' will be observed, see Fig. 
[1] (thick lines). The corresponding excitation probabil- 
ity is (see the Appendix A for explicit expressions of the 
matrix elements in Eq. (fTU]) ) 
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4a2 cos — + n^ cos(r2'T) 



4f A^ + fi'2 + fl^ cos —) sin^ — cos $ 
V 2 / 4 . 



},(11) 



where $ is a combination of the three individual laser 
phases, $ = 0i — 202 + (f>3, and ft' = y/i}'^ + A^, see Fig. 
12 (solid line). 

Remark 1: Unlike the Ramsey configuration, this pat- 
tern is independent of T and thus of the intermediate dis- 
tance L between the pulses since both interfering paths 
spend the same amount of time in the upper level and 
there is no accumulation of a phase difference. 

Remark 2: dPgf/d^ ex sin<I> so that there is a mini- 
mum (zero) at $ = independently of all other param- 
eters. In particular, this is true (within this approxima- 
tion) regardless of the velocity and the precision with 
which the tt/2 and n pulses are implemented. 

Remark 3: The detuning may affect the visibility of 
the fringes but does not shift the fringe pattern. Near 
resonance, A <C ri, and considering perfect pulse areas, 
(fir = tt), Eq. pT|) (which does not depend on these 
conditions) may be expanded to leading order in A as 



P: 
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ge 
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^"^-'^-T^ 



(12) 



The semiclassical central zero is in summary robust ver- 
sus velocity or detuning variations, and does not require 
strict (7r/2 and tt) conditions on the pulse areas. How- 
ever, if the kinetic energy of the atom is comparable with 
the interaction energy, this semiclassical approach breaks 
down and a full quantum mechanical solution becomes 
necessary, yielding a phase shift of the interference pat- 
tern as we shall see. 



IV. QUANTUM TREATMENT 

The general solution to the stationary Schrodinger 
equation ([4]) away from the laser fields takes the form 

1^.) = (G+e'^-^'- + G_e-''=^-)|5) 

+ (S+e'9^'^" + S_e-'''-")|e), (13) 

where the amplitudes G± and E± have to be determined 
from the boundary and matching conditions. We will 
follow Ref. [2l| to derive the exact quantum result of the 
interference pattern. Let us denote by Rlj {T^j ) the total 
refiection (transmission) amplitudes of an atom entering 
the interferometer from the left in the i channel and an 
outgoing plane wave in the j channel, i,j — g,e. In the 
case of a plane wave incident from the left in the ground 
state, we have, for the leftmost laser-free region I, see 
Fig.m 



1^). 






Jll^^e'^k^^-j 



\9) 



Rge^ 



-iq^x\ 



(14) 



whereas the outgoing wave to the right of field 3, region 
VII, has the form 



„=r' e'^--|5)+T'e' 



•I „iq^x\ 



(15) 



That is, after passing the three laser pulses the atom 
may still be in the ground state, propagating with a 
wavevenumber k^, or in the excited state, propagating 



with a wavenumber qx- In the latter case, the atomic 
transition \g) — + |e) induced by the laser field changes the 
kinetic energy in the effective equation for x-direction. 
For A > the kinetic energy of the excited state com- 
ponent is enhanced by ?iA whereas for A < it is re- 
duced by h/S.. For A smaller then the critical value 
Acr = —hk'^/2m, the excited state component becomes 
evanescent and its transmission probability vanishes (the 
channel becomes closed). Thus, the quantum mechan- 
ical probability to observe the transmitted atom in the 



excited state is zero for A < A^^; otherwise 
pg _ fe I ; ,2 A > A 

Se ~ U I gel iUl ^ ^ l-^cr. 



(16) 



and we shall limit the analysis to this later case. (For a 
study of the evanescent regime see [23 )• The exact form 
of r' follows from the matching conditions between the 
free-space solutions and the dressed state solutions inside 
the fields using the transfer matrix formalism J2ll. l26l |27| . 



A. Excited state probability amplitude 



The solutions (|T3|) in the laser free regions may be given 
in a compact form by a constant 4-dimensional vector 
V — {G+,G-,E+,E^)' (the prime means "transpose") 
with the complex amplitudes. In particular, the scat- 
tering boundary conditions are imposed on the external 
regions / and VII, 



(l,i? , 0,i? )', 



r' 

99' 'ge 



V = fT' r' 0)' 



(17) 
(18) 



These solutions at the external regions are related by a 
combination of transfer matrices [2ll |27[ , see Appendix 
IE as 



Eq. (|2T|) provides an explicit, and easy to calculate ex- 
pression. However, this numerical calculation alone does 
not necessarily provide much physical insight. It is useful 
to relate the transfer matrices to scattering transmission 
and reflection amplitudes for the individual laser regions. 
We denote by (r„)'j and {tn)\j the single-laser reflection 
and transmission amplitude for incidence on the nth laser 
"barrier" from the left in the ith channel and an outgo- 
ing plane wave in the jth channel (as before i,j = g,e), 
and by (?"„)[■ and {tn)ij the corresponding amplitudes 
for right incidence. These scattering amplitudes for the 
laser units are also easy to calculate (exactly, using Eqs. 
(IC6IC8IB10|) . or with approximations, e.g. semiclassi- 
cally [2y| or otherwise), and their moduli are typically 
close to one or zero, so that we may introduce an ex- 
pansion parameter, see below, to discern the dominant 
contributions corresponding to "direct scattering", and 
classify the order of the corrections in terms of the num- 
ber of reflections. One further advantage is that we may 
also classify and distinguish the paths according to the 
number of inter-laser free-motion regions in which the 
atom flies in the excited state (for direct, reflectionless 
paths, this number may be 0, as in A4 of Fig. [H 1, as 
in A2 and A3, or 2, as in Al). We may thus distinguish 
those paths that will finally interfere (e.g., A2 and A3 in 
Fig. [fl from those that wiU not (Al and A4 in Fig. [T|) 
in an atomic wave-packet because of the arm separation 
due to recoil. 

(n) 

The relation between the matrix elements T, ■ and 
the individual scattering amplitudes (?'n)i'w, (^n).'!; is in- 



vertible. 



1. 



t(") — f 
' a — Jn 



Vn)ij, [inji j 



known functions /„ (see Appendix [Cjl . 



2, J ' \ "'''2 J 

with invertible 



Mach-Zehnder terms 



^.=tWt(^^T(3)«.„=T-*^.„, (19) 

where T^"^ is the transfer matrix for the nth laser, i. e., 

T(i) = T(O,Z/2,0i), 

T(2) = T(L + V2,L + 3V2,02), 

T^^) = T{2L + 3l/2,2L + 2l,(t)3), (20) 

and T*°* is the "total" transfer matrix for the whole in- 
terferometer. Solving Eq. (fT9| for Tjg, we find that the 
excited state transmission probability amplitude is given 
by 



rpl 



-r-tot 
' 31 



-r-tot-r-tot 

' 13 '31 



-r-tot-r-tot ' 
'11 '33 



(21) 



B. Two-channel recurrence relations 

Since the transfer matrices for the laser interactions are 
known in terms of the laser parameters, see Appendix iBl 



If the kinetic energy {Tikx)"^ / {2m) is larger than the 
Rabi energy hVi the scattering process will be dominated 
by transmission through all fields and all reflection am- 
plitudes will be small quantities compared to the trans- 
mission ones, i. e., \{tn)il\ ^ l('^n)if I fo^' all lasers 
[n = 1,2,3). Moreover, the second laser is assumed 
to apply very nearly a 7r-pulse which flips the internal 
atomic state, so that 1(^2)^^ I ^ 1- Multiplying all small 
amplitudes by a small expansion parameter 77, the series 
expansion in 77 of T' has the form 



22 



Tl = {Al +Al) + r, {A\ + Al) + ry^ ^ Sf , 



(22) 



where the individual quantum amplitudes Al and Bl are 
given in AppendixlDlin terms of the single laser scattering 
amplitudes (r„)[j and (i„)[j . 

The dominant zeroth order terms are in correspon- 
dence with the A2^^ and A^^ direct paths in the semiclas- 
sical picture, whereas first order corrections correspond 




FIG. 2: P^f (Eq. ([TT|, solid line) and P«<, (Eq. (JSTJ, dashed 
line) as a function of the phase difference $ for a resonant 
condition A = 0. Data: v^ — 1 cm/s, I — lO/^m, mjva = 
3.82 X 10~^® Kg. The Rabi frequency is fixed to satisfy the 
TT-pulse condition, f2 = nVx/l = 27r x 500 Hz. 



to the semiclassical Af^' and ^|^' paths and are small 
since they contain diagonal transmission amplitudes in 
the second laser (Appendix ID 1[) . The quantum reflec- 
tion effects are included in the second order terms, which 
contain the quantum amplitudes B^ for all the 22 possi- 
ble paths leading to a transmitted excited atom including 
two reflections, see Appendix ID 21 

In a wave packet, recoil effects will separate in space 
all these paths leading to transmitted excited atoms. In 
order to compare the quantum interference pattern with 
the one obtained semiclassically in Eq. pT|) . we choose 
those paths interfering with A^ and v4|, i.e., the ones in 
the Mach-Zehnder geometry (along thick lines in Fig. [1]): 
they are B'^ with i = 1,2,3,4. This gives the following 
transmission amplitude. 



i=l 

and the excited state probability 

Al + Al + Y^B! 
see Eq. ([TO]). 



pq ^ Si. 



Let us first write the amplitudes in terms of their moduli 
and phases. 



Al = {ti)Ut^)i,{h)U 



(26) 



l(^i)L(i2)L(i3)' |e-(^ 



l-4'2+<l>3)„i02 



At - (^l)L(^2);,,(i3)L = l(tl)L(^2)^e(i3)L|e-*e^'^ 

where the 6*2,3 transmission phases result from the addi- 
tion of the phases of the individual transmission ampli- 
tudes along the path when all (pn ~ 0. From the semi- 
classical expressions. 
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ASCl 



\A 
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AT^le'^'^'^e''^^'' 



(27) 
(28) 



we have that, in a 7r/2— tt— 7r/2 configuration with A := 0, 
6*1^' — Of^'' — TT, see Appendix A. For the quantum case 



we may also expect 
actual phase as ( 



TT. If we write the 
TT -f (5$, there will be a minimum of 



\Al + Al\' = \Al\^ + \Al\' + 2 cos{^^ + 



Se)\Al\\Al\, (29) 



at $ = (5<i>, 






oc sin (-$ + TT + 5$) = 0. 



(30) 



This is a quantum phase shift which vanishes in the semi- 
classical limit. Imposing the condition Ql/v^ = tt at 
each velocity (i.e., exact semiclassical 7r/2 and tt con- 
ditions), the quantum motion shift is shown in Fig. [2] 
(dashed line). Take note that the minimum of P^g is 
not a zero since the quantum moduli jA'l and |A|| do 
not exactly coincide: these conditions do not really split 
the beam in two equal halves at the external lasers, so 
that \{t,V\ ^ \{t,y\, and \{t,)L\ ^ HhYj. The con- 



sequence is a quantum reduction of visibility. We may 
look for enhanced visibility modifying the laser intensity 
and thus the pulse area away from the former condition, 
i.e., Vtl/vx = TT + t. Figs. [3] and 2] show the moduli \AW 
and 1^3 1 and the phase shift 5^ as a function of the extra 
phase e for fixed laser width and velocity. Note that the 

(23) moduh of \AW and \AW cross each other at a value Eq, so 
that a zero of P^g can indeed be achieved by adjusting 
the Rabi frequency at (tt -|- eo)vx/l- There are however 
two important differences with respect to the semiclassi- 
cal exact 7r/2 — tt — 7r/2 case: (a) the "optimal" value Eq 
depends on the velocity (in the semiclassical case £0 ~ 

(24) for all v^)] (b) even for the optimal eo, there is a quantum 
phase shift, (5$ ^ 0. 



E. Reflections 



D. Direct scattering and quantum shifts 

We shall first work out the direct scattering case in 
which the reflection terms can be neglected. 



P',e^^\Al+Alt 



(25) 



Adding the next order in Eq. ([M]) . and taking into 
account the phase dependence of each of the quantum 
scattering amplitudes (Appendix ID 3|) , one may write 



pq 

ge 



qx_ 



-i^ 



{ai 



Bl 
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Al 
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(31) 




FIG. 3: 1^42! and [^43! as a function of e (the deviation from 
perfect pulse areas, i. e., Ql/v^ = 7r+e) for different velocities: 
The thin lines correspond to u^, — 0.5 cm/s and the thick ones 
to Uj, = 1.0 cm/s {I = 10 /xm in both cases). Note that the two 
moduli cross each other at some value of e. At these values, 






A\ I and maximum visibility will be obtained, see Eq. 
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FIG. 4: Quantum phase difference for non-perfect pulse areas 
for fixed values of v = 1 cm/s and Z = 10 /xm. Note that in the 
semiclassical case this difference is 82^^ — 0^'^' = -k for every 
value of e, i.e, for every pulse length. 



where the tildes are the amplitudes for 0„ — Q^ n — 
1,2,3. P^^ is plotted in Fig. [H for ni/v^ = tt. For a 
velocity region the result cannot be distinguished from 
the calculation with direct paths only, which should be 
dominant for h fc^/2m ^ fiVl. Combining this with the 
TT-pulse condition, the direct scattering approximation 
is valid when k^l ^ 27r, as it is observed in Fig. [5l 
where, for lower values of kxl the direct approximation 
breaks down and quantum reflection effects become rel- 
evant. The effect is a rather chaotic oscillation of the 
shift (the actual structure is even more complex than the 
one shown in the scale of the figure) . There are however 



several reasons why this regime will be difficult to see in 
practice as commented in the final discussion. 



F. Effect of the detuning 

The calculations so far have been made for perfectly 
resonant interactions, i. e., for A = 0, where the detuning 
A contains both the natural detuning Ag and the kinetic 
detuning, see Eqs. (|6l7p . In a wave packet it is not 
possible to fulfill the perfect resonant A = condition 
exactly for all components: Even though one may adjust 
the laser frequency to compensate for the recoil term in 
Eq. ([7]), the momentum spread of the wave-packet in the 
y direction, Afcj,, will lead to a detuning spread from the 
Doppler term, A^) — hkykL/m. In the semiclassical case, 
we have already shown that the detuning can only affect 
the visibility of the interference fringes but will not affect 
their position, see Eq. (fT^ . This is no longer true in a 
fully quantum calculation. Since the position spread in 
the y direction, Ay, cannot be larger than fikLT /2m in 
order to suppress the interference with outer paths, Afcj, 
will also be limited. We may thus estimate the Doppler- 
detuning spread as 



A^ 



hkLi^ky 



2 



±^. (32) 



Numerical simulations with v^ ~ I cm/s and I — 0.1 — 1.0 
m show that a kinetic detuning like this has negligible 
effect in the calculated phase shift, which is quite robust 
against detuning ffuctuations. 



V. DISCUSSION 

In this paper we have explored the low velocity limit of 
atomic interferometry in a simple Mach-Zehnder 7r/2 — 
tt — tt/2 configuration of spatially separated laser fields ig- 
noring further external fields. In particular, we have per- 
formed a fully quantum analysis of incident monochro- 
matic stationary atomic waves by providing explicit ex- 
pressions for transmission probabilities from which the 
physically relevant paths and contributions in terms of 
transmission and reflection amplitudes for the individual 
laser fields may be extracted. 

For laser fields separated in space, the ideal 7r/2 — tt — 
tt/2 conditions leading to perfect splitting, perfect reflec- 
tion, and interferometer phase given exclusively by the 
laser field phases cannot be reached in a fully quantum 
scenario, even for a fixed incoming velocity. The conse- 
quence is a quantum-motion phase shift at low atomic 
velocities related to the phases of the transmission am- 
plitudes. One may optimize the fringe visibility by de- 
viating the Rabi frequency from the semiclassical value, 
but a phase shift remains which, in addition, depends on 
the incident velocity. This quantum-motion shift is quite 
insensitive to the detuning to be found in wave packet 
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FIG. 5: Phase shift of the interference pattern as a function 
oi kxl. The lenght / is kept constant, while the velocity (and 
thus the Rabi frequency) are changed in order to maintain 
the 7r/2-pulse condition. It should be clear from Eq. (|lip 
that in the semiclassical regime, there is no phase shift (solid 
line). The dashed line corresponds to the direct scattering 
approximation, where quantum reflections are neglected. At 
low velocities, quantum reflections become relevant and direct 
approximation breaks down (dotted line). Resonant A = 
pulses have been considered for the calculations, but numeri- 
cal simulations show the robustness of the phase shift against 
detuning fluctuations. Changes in the detuning of the order 
±2vx/L are indistinguishable in the scale of the figure, see 
Sec. ITVFl 



components but shows wild oscillations when the veloc- 
ities are so low that paths with reflections at the fields 
become significant. 

All the above has been done for square laser profiles in 
the longitudinal direction, with two-channel recurrence 
relations which are by construction well adapted to gen- 
eralizations for more realistic laser intensity profiles. We 
have also considered flat laser sheets ignoring the cur- 
vature of the field. For direct paths (transmitted in all 
lasers) this is a good approximation, whereas paths with 
refiection, having longer flights and more collisions with 
the laser fields, will be more affected by curvature effects, 
which, together with other averaging effects (because of 
their extreme sensitivity to tiny velocity variations) will 
surely cancel their contribution to the shift. 
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APPENDIX A: AMPLITUDES FOR ATOM AT 
REST 

These are the amplitudes needed for calculating matrix 
elements in Eq. (jlOp with the Hamiltonian ([5]). Here the 
atom is illuminated during a time t. 

{9\e-'"-"%) = e*^*/2[cos(r!'t/2)-*Asin(fi't/2)], 



^' 



A 



(g|g-«ff„tA|g) ^ e''^''^[cos{Q.'t/2) + i^s\n{Sl't/2)] 



^' 



Vt 



{e\e~'"-''%) = -ie*^*/2e-**"— sin(r!'V2), 



(^|g-^H„t/'i|e) = -ie*^*/2e^'^"^sin(17't/2), (Al) 



0/ 



where J7' = vfP + A^. They are also useful to obtain 
semiclassical approximations of matrix elements of trans- 
mission amplitudes i„. 



APPENDIX B: TRANSFER MATRICES 

Consider the regions a — I, II, III in Fig. [I] separated 
by xi and X2- The general solution to the stationary 
Schrodinger equation (j4]) of the effective Hamiltonian ([S]) 
reads 



\'ijj{x))a = gaix)\g) + ea{x)\e) 



(Bl) 



We want to find these solutions for a — I, II, III and 
match them at the boundaries. 



1. Solution outside and inside the fields 

The solutions at the laser-free regions {a = /, ///) are 
given by 



|V(x))„ = (flae^'^-" + 6ae-^'=-^) |5) 



+ (^ 



d„e 



(B2) 



where hkx is the initial momentum of the atom in the 
longitudinal x direction and q^ — k'^ + 2mA/'h. Inside the 
laser fields the (unnormalized) dressed state basis which 
diagonalizes the interaction part of the Hamiltonian is 
given by |A±) = \g) + 2\±e-''^n-^\e), where X± = (-A± 
Q!)/2 are the dressed energies. The solution inside the 
interaction region {a = II) will be given in terms of these 
dressed states and dressed energies, 



+ (-^ 



Cue 



ik-2 



iiie 






)|A_), (B3) 



with wavenumbers fcj. = /c^ — 2m\±/fi. The solution in 
each zone can be then given by a set of 4 unknown com- 
plex amplitudes, collected in a constant complex vector 
Va = {aa,ba,Ca,da)' , whcrc the prime means "trans- 
pose" . 



2. Matching Conditions: one laser 

The wave functions and their derivatives with respect 
to X may be written in the foUowing way in each of the 
zones. Outside the interaction region (a = /,///), 



(B4) 




(B5) 



where the dot represents derivative with respect to x. 
The M matrices are exphcitly given by 



Mo(x) = 
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^ik+x 
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n 



-e ^e 

ik^ 



p — ik^x 






2A— ^^icf) ^ik^x 
n ^ ^ 

ik-X 



2A- ^ — i(f)^ — ik-x 
n ^ ^ , 



\ 



ik+e +^ ~ik+e +'^ ik^e 

\ik+^^e-''t'e^''+=' -ik+^±e-^'t'e-'''+'' ik^^e-"t'e^''-'= -ik^^^e-^'t'e-'''-'' J 
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(B6) 



With this notation, the matching conditions at a; = xi 
and X — X2 can be written as 



Mo{xi)Vj = Mb(xi,0i)v„, 



(B7) 
(B8) 



Ehminating v ^^ from the system above, we end up with 
a transfer matrix T(xi, X2, (t>i) which connects the aniph- 
tudes of both sides, 

Vj =Tixi,X2,(t>i)Vjjj, (B9) 

defined by 

T{xi,X2,(l)) = Mo{xi)-^Mi,{xi,(j))Mb{x2,(l))-^Mo{x2). 

(BIO) 



4. Multiple laser fields 

Clearly we may repeat step by step the operations 
above for the second and third laser. The results are for- 
mally the same, except for the substitution of the match- 
ing points and the laser phase. We may then write 



Vj = T{xi,X2,4'l)Vjjj, 
'"ill = '^{X3,X4,(P2)V^,, 



(Bll) 

(B12) 
(B13) 



and relate the waves on the extremes by 
Vj ^T{xi,X2,(f>i)T{x3,X4,(j)2)T{x5,xe,(f)3)v^jj. (B14) 



3. Phase dependence 

The explicit dependence of the (one laser) transfer ma- 
trix on the laser phase (we drop the laser index n) is 
as follows 



T(xi,X2,( 



/ tn ti2 e^*ti3 e^-^tu^ 

721 t22 6^*723 e*"^t24 

Ve-^*t4l e'''f'T42 t43 744 / 



where the tildes represent the phase-free form of the am- 
plitudes, i. e., Tij = Tij{(j) = 0). 



APPENDIX C: RECURRENCE RELATIONS 

Consider, for the nth laser located between Xi and x/, 
the following "elementary" scattering boundary condi- 
tions corresponding to incidence of a wave in one channel 
from left or right: 

• Left incoming, ground state: 



Vnjgg 





_t(") 



T^"'{xi,Xf) 





V'njge 

V J 



(Cl) 



Left incoming, excited state: 
/ \ ( {U 



1 



= T(")(a:.,x/) 



V 





V'njee 





• Right incoming, ground state: 



(tr. 







'gg 



Right incoming, excited state: 
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{rnVge 







(C2) 



(C3) 



(C4) 



Thus, for each laser we have a system of 16 equations 
which can be solved to give the transfer matrix T^"^ el- 
ements as a function of the single field scattering ampli- 
tudes {ra)lj , or the other way around, the single field 
scattering amplitudes in terms of the transfer matrix ele- 
ments. Combined with Eq. (jBlOp . this provides explicit, 
exact expressions for the scattering amplitudes. 



with the common denominator / defined by 
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2. (r„) -'J and (tn)''J as a function of T[ 



r{") 



as a function of (r„)/^ and (f„ 



We have dropped the n index of the laser for simplicity. 



We have dropped the n index of the laser for simplicity. 
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''ge 



— T13T22T31 + T12T23T31 + T13T21T32 — T11T23T32 — T12T21T33 + T11T22T33 

F 
T13T32T41 — T12T33T41 — T13T31T42 + T11T33T42 + T12T31T43 — T11T32T43 

F 
T14T23T31 — T13T24T31 — T14T21T33 + T11T24T33 + T13T21T34 — T11T23T34 

F 

~Ti4T33T4i + T13T34T41 + T14T31T43 — T11T34T43 — T13T31T44 + T11T33T44 

F 



(C8) 



where the common denominator F is given by 



second laser flips the atomic state are explicitly given by 



F — T13T31 — TllT; 



33 



(C9) 



APPENDIX D: EXPLICIT EXPRESSIONS OF 
THE QUANTUM SCATTERING AMPLITUDES 



We give here explicit expressions for the amplitudes in 
Eq. dMD- 



1. 4 direct paths 



There are four possible paths leading to an excited 
atom with no reflection. These are the corresponding 
amplitudes, 
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2. 22 paths with two reflections 



3. Phase dependence of the scattering amplitudes 



The quantum amplitudes Bf for all the 22 possible 
paths leading to a transmitted excited atom including 
two reflections, provided that the perfect 7r-pulse at the 



The dependence of each of the path amplitudes on the 
laser phases are easily obtained from the 2-channel re- 
currence relations and the transfer matrix formalism. If 
the phase-free amplitudes (for all 0„ = 0) are denoted by 
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tildes, we have 
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